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ABSTRACT 

A  ganaral  nodal  for  Cha  daslgn  of  float  aslntensnce  systems  Involving 
several  craws  working  In  aaquanca  la  davalopad.  Raphaels  la  placad  on 
optlalzlng  tha  system  for  schadulad  parlodlc  aalneananca  activities*  Results 
ara  glvan  for  tha  easa  where  craw  service  tlaas  are  exponentially  distri¬ 
buted.  When  naan  service  tines  for  all  crews  ara  equal,  an  optlaal  systea 
configuration  Is  deterninad  as  a  function  of  the  nueber  of  craws,  cost  of 
craw  activities  and  cost  of  a  fleet  unit's  Idle  tlaa. 
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A  FLEET  MAINTENANCE  SYSTEM  DESIGN  MODEL 
INTRODUCTION 

Fleet  maintenance  lyitna  typically  arc  designed  to  service  fleet  units, 
such  as  trucks,  airplanes,  cars,  etc*,  on  a  scheduled  maintenance  program  as 
well  as  emergency  repair  basis*  Basic  system  design  variables  include  the 
division  of  work  between  service  crews,  level  of  service  provided  by  such 
crews  and  a  limit  on  the  number  of  fleet  units  idle  at  one  time* 

A  general  model  for  the  design  of  such  systsms  is  developed  and  described 
in  this  article*  Since  in  a  given  fleet  environment  the  coot  end  effect  of 
emergency  repairs  tend  to  be  independent  of  system  design,  emphssls  is  placed 
on  optimising  the  system  for  scheduled  maintenance  activities* 

Typical  fleet  melntenance  systems  ara  made  up  of  multiple  crews  and 
skills  which  must  be  scheduled  in  a  fixed  sequence  to  complete  specific 
maintenance  tasks  on  each  unit.  Task  times  usually  vary  considerably  from  one 
unit  to  the  next.  Studying  such  systems  is  import ent  because  of  the  large 
number  of  organisations  confronted  with  this  problem  and  the  magnitude  of 
investment  and  expense  consumed  by  such  systems* 

The  pertlculsr  model  developed  here  considers  s  fleet  of  erbltrsry  slss 
and  a  service  fecillty  consisting  of  a  finite  number  of  crews  working  In  a 
fixed  sequence  on  each  fleet  unit  removed  from  operation*  A  limit  is  placed 
on  the  number  of  units  which  can  be  removed  from  the  operating  fleet  at  any 
one  time*  Once  the  limit  la  reached,  the  next  scheduled  fleet  unit  must  wait 
until  another  unit  is  releesed  from  the  maintenance  system  before  it  can 
enter.  In  its  general  form,  system  coats  ara  a  function  of  scheduled 
maintenance,  emergency  repairs  and  opportunity  or  capital  costs  associated 
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with  14U  oparatlag  wait*.  Capital  eoata  and  oparatlag  expanses  associated 
with  heaolag  tha  ayataa  ara  not  considered*  In  tha  gaaaral  nodal,  craw 
aarviea  tlwaa  ara  eoaaldarad  to  ba  ladapaadaat  randow  variables  with  a  gaaaral 
distribution*  la waver,  tha  analysis  hara  la  raatrlccad  to  tha  caaa  whara 
aarwiaa  tlaaa  ara  aapoaaotlally  diatribuead.  When  aaaa  service  tlaaa  for  all 
erowa  ara  aqwal,  aa  aptlaal  ayataa  configuration  la  dotoralnod  aa  a  faaetloa 
of  tha  aoahar  of  aiawo,  tha  ooot  of  craw  activities,  aad  tha  coat  of  a  float 
wait'a  idla  tlaa*  inch  eoafiguratlona  ara  daacrlbad  la  taraa  of  tha  Halt  on 
tha  anabar  of  waita  which  caa  ba  taasvad  froa  tha  oparatlag  float  at  aa y  oao 
elaa  aad  tha  level  of  service  rat a  necessary  for  aaeh  craw* 

GENERAL  COST  MODEL 

Coaaldar  tha  following  aodal* 

N 

Z  -  B[P(R(j»,y))J  ♦  Y[L(^,y)]  ♦  l  CCO 

kvl 

whara 

Z  *  expected  total  systan  coat  par  unit  tlaa  (In  steady  atata) 

K*)  ■  axpactod  coat  aaaoclatad  with  aaargancy  rap air a  par  unit  tlaa 

t(*)  -  axpactod  coat  aaaoclatad  with  float  units  ldla  during  thair 

raaoval  froa  tha  oparatlng  float  par  unit  tlaa 

C(p^)  •  axpactod  coot  aaaoclatad  with  aalntananca  craw  k  par  unit  tlaa 
aad 

F(«)  •  axpactod  nuabar  of  units  r squiring  aaargancy  rapalra  par  unit  tlaa 

t( • )  *  axpactod  production  rata  (rata  of  coaplatlng  aalntananca  of  units) 

of  aalntananca  ayataa 

*  (S(6jJ)~*,  tha  service  rata  for  craw  k  • 
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where 


S{6k)  -  the  expected  carries  Clac  for  crew  k 

ji  •  the  rector  whose  components  ere  the  ^(k  -  1,2,3,...,  N) 

H  ■  the  number  of  erewc  In  systaa 

Y  -  nl— i  nuaber  of  fleet  units  reaoved  froa  operating  fleet  at  any 
point  In  tiae  which  ere  available  to  maintenance  crews 

L<«)  •  expected  nuaber  of  fleet  waits  In  the  aaintananee  system. 

Components  of  cost  considered  here  relate  to  emergency  repair, 
opportunity  costs  of  Idle  units  and  crew  costs  associated  with  scheduled 
maintenance  activities.  It  is  recognised  that  actions  can  be  taken  and,  where 
possible,  should  be  taken,  to  reduce  •<•);  however,  for  our  purposes  KO 
Is  considered  to  be  constant  and  so  does  not  Influence  the  design  of  the 
aaintananee  systaa  required  to  service  scheduled  aaintenaaee  activities • 

SEQUENTIAL  CUV  MODEL 

The  problea  now  becomes  that  of  solving  for  the  values  of  |ip  y^,  ... , 
nB  and  y  that  alnialxa  total  systaa  cose,  where  the  chosen  value  of  each 
pfc  then  would  be  (approxlaatsly)  realised  by  adjusting  the  slse,  composition, 
end  node  of  operation  of  crew  k  appropriately.  For  those  systeas  which  ere 
kept  operating  at  as  does  to  capacity  aa  possible,  given  the  values  of  these 
design  variables,  the  Sequential  Crew  Model  (1]  can  be  used  to  provide  systaa 
production  rates  as  a  function  of  the  nuaber  of  crews,  crew  service  rates,  and 
the  Halt  on  the  nuaber  of  fleet  units  available  to  the  aaintananee  systaa. 


Th«  assumptions  node  by  chls  aod«l  art  tha  following 


(1)  Tha  syseaa  of  N  crews  la  available  continuously. 

(2)  Thaaa  H  craws  always  parfora  chair  tasks  on  aaeh  unit  in  eha  saaa 
fixed  ordar,  wharo  eha  naxc  craw  begins  (if  noc  already  occupied  on 
anothar  unit)  aa  soon  as  eha  procadlng  craw  anda. 

(3)  For  crew  k,  eha  service  elass  for  eha  raspacelva  units  are 
independone  and  identically  dlacrlbuead  according  eo  aa  exponential 
distribution  wleh  paraaacer  • 

(4)  Ho  wore  chan  y  unlea  can  be  "in  eha  syseee"  (already  begun  by  at 
lease  eha  first  crow  but  not  yet  finished  by  eha  lose  craw)  at  one 
else. 

(5)  Escape  whan  eha  syseaa  la  at  this  upper  llaie,  eha  firae  craw  always 
scares  work  on  another  unie  aa  soon  as  it  coaplaees  its  work  on  its 
currant  unie. 

Tha  following  resules  than  are  obtained  frow  eha  Sequential  Craw  Modal. 

t<£.Y>  -  p  l  f. 

1  lag  1 


U)»,y)  -  I 
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where 


-  stationary  probability  of  being  in  seaca  1  of  eha 
underlying  Markov  chain. 

■  nuabar  of  unles  in  eha  syseaa  in  seaca  i  . 

S  •  eha  see  of  P^  wleh  eha  flrse  craw  busy. 


A  more  complete  description  of  this  model  and  let  chtrtcCtrLaelc  bthtvloc 
eta  bt  found  In  [2]. 

le  ctn  bt  readily  obatrvad  that  Chit  nodal  can  bt  fornulaead  at  a 
coneinuoua  tint  paraatttr  Markov  chain;  hence,  eha  can  ba  obealnad  from 
cha  solution  of  a  ate  of  llnaar  aquations.  Slnca  aost  practical  systaas 
raqulra  larga  sacs  of  aquations,  eha  raadar  aay  wish  to  rsftr  to  (5]  and  [6] 
for  procaduras  for  simplifying  solution  asehods. 

Moans  art  now  available  to  datoralna  an  optiaal  design  by  calculating 
lt(*)  and  L(0  and  using  an  appropriate  numerical  starch  routine.  An 
example  using  a  direct  nuaerleal  approach  Is  given  in  Appendix  A  for 
Illustrative  purposes. 

BEHAVIOK  OF  BALANCED  SYSTEMS 

A  very  useful  and  frequently  obaarvad  case  of  this  system  occurs  when 
expected  service  times  of  all  creua  are  equal,  l.a.,  •  |i  for  all  k  .  An 

examination  of  the  transition  Intensity  matrix  for  this  case  reveals  that  It 
Is  doubly  stochastic.  Hanca , 

P1  “  m  £or  *11  1  » 

where  m  is  the  number  of  seates  In  the  Markov  (Stain  and  Is  the 
stationary  probability  of  eha  system  being  In  state  1  (sea  [4]). 

This  fortunate  relationship  makes  It  possible  to  calculate  P^  by 
counting  the  number  of  staeas  In  eha  Markov  Chain  representing  a  particular 
syseem.  Such  systems  and  respective  Markov  Chains  are  defined  by  a  given 
number  of  crews,  crew  service  rates,  and  a  Hale  on  the  number  of  units 
allowed  In  the  system. 


OPTIMIZING  BALANCED  SYSTEMS 


a 


I  jH'i 


Derivations  of  R(*)  and  L(*)  for  balanced  systems  are  given  in 


Appendix  B.  For  such  systems 


«•>  •  “  •  (ffe) 


uo  .  y  -  Jcix=12 

1  y  NCy+N-1)  * 


where  p  is  Che  common  value  of  the  pk  (k  -  1,2,...,N). 


aw 

Now  assume  that  C(ji)  -  J  C(pk)  is  an  increaaing  linear  function 


of  p  , 


C(ji)  -  cp  where  c  >  0  , 


and  that  Y(L(*))  Is  an  increasing  linear  function  of  L( •)  , 


Y(L(»))  ■  ML(»)  where  M  >  0  . 


Further  assume  that,  for  a  fixed  integer  N  >  I,  the  objective  is  to  select  a 
(positive  Integer)  value  of  y  so  as  to 


minimise  Z  •  cp  +  ML(j*,  y) 


subject  to 

R(li.t  r)  -  B*  . 


Thus,  R*  is  the  press cablished  mean  race  at  which  routine  maintenance  of 
fleet  units  would  be  scheduled*  (Note  ehat  R*  is  the  mean  rate  rather  than 
the  actual  rate  over  any  short  period  of  time,  since  the  actual  rata  at  which 
units  eneer  the  system  is  p  when  the  number  already  in  the  system  is  less 
than  y,  whereas  this  rate  becomes  sero  while  this  number  equals  y.)  Given 


Vk 


the  mnb«r  of  units  la  Cho  fleet,  tho  value  ehoaoa  for  R*  normally  would  bo 
baaod  oa  tho  doolrod  aalacoaaaco  eyelo  (clao  botwooa  aehodulod  malntonaneo 
events)  for  tho  Individual  ualca. 


For  any  given  cholco  of  y ,  satisfying  tho  eoaatralnt  of  tho  problem 
onouata  to  chooalag  tho  value  of  p  that  equatea  R(ji,y)  with  R*,  namely. 


(*) 


Therefore,  Who  problem  (?)  can  be  reatated  aa  aelectlng  a  (poaltlve  integer) 
value  of  y  ao  aa  to 

(?*)  mlalalae  Z  -  cR*  +  H[y  -  . 


To  begin  aolvlag  thla  problem  (?'),  .the  Integer  restriction  shall  be 
dropped  temporarily  In  order  to  treat  y  aa  a  contlnuoua  variable.  Call  thla 
relaxed  problem  (P).  After  conelderable  algebra,  the  first  two  derivatives 

can  be  expressed  as 


cR*(K-l)  .  „ 
2  *  “ 


M  (H-l)M 
«T  11  2  » 

(y+H-l) 


-  2(H-i)  - - — x  1  . 

6y  y  <y•«l-l)', 


The  common  approach  now  to  finding  a  global  minimum  would  be  to  set  the  first 
derivative  to  aero  and  solve.  However,  since  Z  la  not  a  convex  function 
of  y  la  general  (note  that  the  second  derivative  becomes  negative  for 
sufficiently  large  y  If  H  la  large  relative  to  cR*),  further  analysis  la 
needed  to  show  that  this  necessary  condition  for  optimality  also  is 
sufficient.  This  la  done  by  proving  the  following  theorem. 


-'*•»*  ^  ■  ♦'JLnrv  ftA  H»v  .<s«»  J^e*. 


Theorem  1:  For  problem  (P) ,  chare  exists  e  unique  positive  y  •  y*  such 


<  0,  if  0  <  v  <  y* 


-  0,  if  y  -  y* 


>  0,  if  y  >  y*  , 


so  y  *  Y*  achieves  the  global  minimum. 


Proof:  Since  N  2,  it  is  clear  that 


Also  note  chat 


jr”  <  Of  ““!r  >  0  for  sufficiently  small  y  >  0  . 

°y  6y 


14.  6*  .  M  M  S  ft 

M  H  >  0  • 


Since  is  a  continuous  function,  it  thereby  assumes  the  value  zero  at 

some  point  by  the  Intermediate  Value  Theorem.  The  proof  now  reduces  to 


O  6 

since  this  implies  that  ^  can  never  return  to  zero  (or  less)  as  y 

6  2  52- 

increases  beyond  the  first  zero  point  of  because  is  a 
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continuous  function  and  so  must  be  strictly  positive  when  is  strictly 
positive  but  sufficiently  close  to  zero. 


showing  that 


SSSETSMOa 


i 
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More  generally,  note  thee 


cR»(N-l)  ,  (N-l)M  <Q 


(y+N-1) 


2  ^cR*(N-l)  _  (N-l)M  j  > 


(y+n-1)4 


2cR»(H-l)  _  (N-l)M 
Y3  (Y+N-1)3 


>  0 


62Z 


6y' 


>  0  , 


which  completes  the  proof. 


Corollary:  If  y*  is  Integer,  it  is  the  optimal  solution  for  (P*)  or 
(P).  If  y*  Is  not  integer,  then  this  optimal  solution  is  obtained  by 
rounding  y*  up  or  down  to  an  Integer  and  selecting  the  one  with  the  smaller 
value  of  Z  for  (P*)  . 


Unfortunately,  solving  directly  for  y*  Is  not  easy  because  it  Involves 
solving  a  quart ic  equation.  However,  the  following  theorem  gives  an  upper 
bound  on  y*  that  is  relatively  tight  when  M  is  small  relative  to  cR*  . 
This  bound  then  provides  the  basis  for  an  efficient  algorithm  for  solving 
(P')  or  (P)  ,  as  summarized  In  the  corollary  to  Theorem  2. 


Theorem  2:  The  value  y  •  Y*  Identified  in  Theorem  1  satisfies  the 
inequality, 

y*  <  a  , 

where 


Proof?  Sine* 

§£  >  (N-l)lf  -  for  *11  y  >  0  , 

it  laaedlately  follow*  froa  Theorem  1  that  y*  1*  bounded  above  by  th* 
positive  root  of  th*  aquation, 


which  yields  th*  d«cir*d  result. 

Corollary:  Lee  A*  b*  th*  l*a*e  int*g*r  greater  than  or  equal  to  A,  and 
let  Z(n)  be  the  value  of  Z  for  (P* )  when  y  •  n  •  An  optiaal  solution 
for  (P')  or  (P)  then  is  identified  by  the  following  aiaple  algorithm: 

(0)  Set  n  ■  A*  • 

(1)  If  n  •  1,  then  y  •  1  la  opeiaal  so  stop. 

(2)  Calculate  Z(n)  and  Z(n-l)  . 

(3)  If  Z(n)  <,  Z(n-l),  then  y  •  n  is  optiaal  so  stop.  Otherwise, 

reset  n  ■"  o-l  and  return  to  Step  1  for  another  iteration. 

Once  the  optiaal  value  of  y  has  been  obeained,  the  corresponding  optiaal 
value  of  n  la  obtained  by  plugging  this  value  of  y  into  equation  (*). 

A  numerical  exaaple  illustrating  this  algorithm  is  presented  in 


Appendix  A, 


It  has  been  uiuMd  so  far  that  tha  nuabar  of  crews  N  has  baan  fixed  in 
advance.  If  M  actually  is  a  dasign  variable  In  addition  to  y,  than  an 
overall  optiaal  solution  for  both  N  and  y  can  ba  obtainad  by  rapaatlng  tha 
above  algorithm  for  aach  of  tha  savaral  possible  valuas  of  N,  and  than 
choosing  tha  ona  (along  with  Its  optiaal  y)  that  ylalds  tha  ovarall  alnlaua 
for  Z* 

Conclusions 

The  usa  of  tha  sequential  craw  nodal  wakes  It  possible  to  analyse  and 
evaluate  systems  performing  maintenance  on  fleet  units*  When  such  saquantial 
craw  systems  are  balanced,  than  It  Is  rslatlvaly  straightforward  to  calculate 
optiaal  values  for  tha  limit  on  tha  nuabar  of  units  allowed  In  tha  aalntananee 
systaa  at  any  ona  tlaa  and  optiaal  craw  service  rates* 
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APPENDIX  A 


Kxsapls  Illustrating  Solution  Proeoduro 
for  Optimising  Balanced  Systeas 

Consider  a  float  of  200  units,  each  requiring  progressed  naintenance 
every  ala  aonths.  This  establishes  R*  at  33y  units  to  be  scheduled  per 
aonth  on  the  average.  As suae  it  has  been  determined  that  the  coot  for  keeping 
a  fleet  unit  Idle  would  be  $600  per  aonth,  and  experience  or  work  standards 
Indicate  that  total  crew  costs  would  be  $10,000  per  aonth  If  the  service  rate 
for  each  crew  were  33j  units  per  aonth  and  each  crew  could  work 
Independently  full  tiae  and  without  intarference  froa  other  crewe.  Also 
consider  chat  the  skills  involved  necessitate  three  crewe  working  in  sequence. 
Thus,  using  notation  defined  previously  and  one  aonth  as  the  unit  of 

tiaa, 

N  -  3,  R*  -  33y,  c  -  $300,  M  -  $600  . 

Applying  Theorea  2, 


so  that  A*  *  8.  Row  applying  the  algorltha  outlined  In  the  corollary  to 
Theorea  2  yields  the  following  results: 

Iteration  It  Z<8)  •  $16,180  sod  Z(7)  -  $16,124. 

Since  Z(8)  >  Z(7),  do  another  Iteration. 

Iteration  2:  Z(6)  -  $16,183. 

Since  Z(7)  <  Z(6),  stop. 

T  •  7  Is  opclaal. 


rlcal  calculation  of  these  Z(r)  la  sunwarlsed  below. 


Hence,  eho  optlaal  design  allows  a  asrlaua  of  7  float  units  Idla  at  any 
ooa  tlao  (y  •  7)  and  provides  suf flelant  aanpowar,  ate. ,  la  aach  crow  to 
achieve  a  craw  service  rata  of  42y  units  par  aoath.  Such  a  ays  taw  will 
process  an  average  of  33x>  fleet  units  per  aoath  at  the  lowest  expected  total 


APPENDIX  B 


Derivation  of  R(»)  and  L(*> 
for  Balanced  Systems 

To  detsrnlns  the  expectsd  system  aorvleo  rata,  !(•),  for  a  balaocad  ijitu, 
one  mot  first  determine  *,  tha  total  number  of  states  in  the  Markov  chain. 
These  states  Indicate  hov  many  of  the  y  units  are  with  each  of  the 
respective  craws,  where  any  maker  larger  then  one  for  the  first  crew 
indicates  that  these  extra  units  are  not  pet  In  the  system  because  the  first 
crew  hasn't  started  on  then  yet.  Therefore,  n  is  equal  to  the  number  of  ways 
Y  Indistinguishable  objects  esn  be  partitioned  among  N  cella,  which 
elementary  combinatorial  analysis  (e.g.,  see  [3,  p.  38])  gives  as 

Y  +  *  -  1 

(Note  that  having  more  than  one  object  at  the  first  cell  corresponds  to  having 
those  units  not  yet  in  tha  system  because  the  first  crew  hasn't  started  work 
on  them  yet.)  The  number  of  states  in  which  the  first  crew  is  busy  is  equal 
to  tha  number  of  ways  (y  -  1)  indistinguishable  objects  can  be  partitioned 
among  N  cells  (since  one  additional  object  has  been  preasalgned  to  the  first 


cell),  or 


Since  the  Markov  chain  la  doubly  stochastic,  the  expected  ijicm  service  race 
chon  eon  bo  ealculaced  as 

«•>  - » £■)  - » • 

The  expected  masbor  of  fleoc  unlcs  in  che  nalncenanco  system,  L (•),  eon 
bo  decernlned  In  che  following  manner.  Trivally,  chia  oxpocCod  nuabor  would 
bo  y  if  Y  units  ware  In  che  system  ac  all  clnoa*  However*  Chore  are  Claes 
when  "empty  spaces'*  (loss  chan  r  units)  exist  In  the  sysceo*  froa  che 
cyclic  queues  equivalent  of  Che  sequential  crew  aodel  (see  [7,  pp.  174-177]) , 
it  can  bo  shown  Chat  Che  tocal  nuabor  of  eopey  spaces  which  occur  in  all  the 
states  of  che  Markov  chain  Is  equal  to  the  nuaber  of  ways  H  lndlsclngulsh- 
able  objects  can  be  partitioned  aaong  (y  -  1)  colls,  or 

H  +  V  -  2 


Hence, 


N 


t-.r 


r  >&  »-»<i  rvv  -  v  > -i-*i '"V  *fc" 
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17.  CMSTMauTIOM  STATCMCNT  (a!  *• 


(  Queueing  Theory 
Rellablllty  Theory*. 

v  J 


Maintenance _  _ 

~z>  Preventive  Maintenance 


rr."'^  i  f1;?"  .i  t  Tnasr 


A  general  «dql  for  the  design  of  fleet  maintenance  systeas  involving 
several  cress  wrking  In  sequence  is  developed.  Raphaels  is  placed  on 
optlalslng  the  systei  for  scheduled  periodic  aalntanance  activities. 
Results  are  given  for  age  case  where  crew  service  tiaes  are  exponentially 
distributed.  When  asan  Service  tiaes  for  all  crews  are  equal,  an  optlaal 
systaa  configuration  is  determined  as  a  function  of  the  nuaber  of  crews, 
cost  of  crew  activities  and  boat  of  a Jleat  unit's  idle  tiae. 
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